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Abstract 

We study the leading order finite size correction (Liischer's /i-term) associated to moving one- 
particle states, arbitrary scattering states and finite volume form factors in 1 + 1 dimensional 
integrable models. Our method is based on the idea that the /i-term is intimately connected to 
the inner structure of the particles, ie. their composition under the bootstrap program. We use an 
appropriate analytic continuation of the Bethe-Yang equations to quantize bound states in finite 
volume and obtain the leading /i-term (associated to symmetric particle fusions) by calculating the 
deviations from the predictions of the ordinary Bethe-Yang quantization. Our results are compared 
to numerical data of the Es scattering theory obtained by truncated fermionic space approach. As 
a by-product it is shown that the bound state quantization does not only yield the correct /i-term, 
but also provides the sum over a subset of higher order corrections as well. 

1 Introduction 

The knowledge of the properties of finite volume QFT is of central importance in at least two 
ways. On one hand, numerical approaches to QFT necessarily deal with a finite volume box and 
in order to interpret the results correctly a reliable theoretical control of finite size corrections is 
needed. On the other hand, working in finite volume is not necessarily a disadvantage. On the 
contrary, the volume dependence of the spectrum can be exploited to obtain (infinite volume) 
physical quantities like the elastic scattering phase shifts |TJ [2] or resonance widths [21 0] . 

Finite size mass correction were first derived by Liischer pp. The 1 + 1 dimensional formulas 
relevant to integrable models together with a generalized F-term formula for moving particles were 
obtained in |6j. Finite size corrections have recently become important in the context of AdS/CFT 
correspondence as well. The generalized /i-term and F-term formulas for moving particles with 
arbitrary dispersion relation were derived in [8]. 

Besides the volume dependence of the spectrum itself, finite volume form factors (matrix el- 
ements of local operators) also represent a central object in finite volume QFT. Apart from the 
obvious relevance to lattice QFT they are also important in 1 + 1 dimensional models, where they 
can be used to construct a systematic low-temperature expansion of correlation functions at finite 
temperature [9]. The connection to infinite volume form factors is given by a simple (though non- 
trivial) proportionality factor [TO], [TTJ [12] which is exact to all orders in the inverse of the volume. 
There are however finite size corrections that decay exponentially with the volume and they play 
a crucial role whenever the numerical simulation is limited to small volumes. Moreover, they can 
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(a) Fusion angles 



(b) The mass triangle 



(c) A triple bound state (see 2.2.3. 
for explanation) 



Figure 2.1: Pictorial representation of particle fusions. 



produce huge deviations even in relatively large volumes provided that the exponent is small. This 
work was partly motivated by such an example which can be found in [12] (section 4.1.2). 

In this work we present a method to obtain the leading /i-term associated to arbitrary multi- 
particle energy levels and finite volume form factors in 1 + 1 dimensional integrable models. Our 
approach is based on the idea that the /i-term is associated to the „inner structure" of the particles, 
ie. their composition under the bootstrap program. It is supposed that the leading /i-term is caused 
by a symmetric particle fusion A a A a — > A c . The results can also be applied in nonintegrable models 
for states below the first inelastic threshold. 

The outline of the paper is as follows. 

We begin our analysis in sec. 2 by giving a new interpretation of Liischer's /i-term and ex- 
tending it to describe moving particles. The scaling Ising model serves as a testing ground for our 
calculations: the analytic predictions are compared to numerical data obtained by the Truncated 
Conformal Space Approach (TCSA) developed by Yurov and Zamolodchikov [13]. The extension 
of our results to arbitrary multi-particle scattering states is presented in section 3. 

Section 4 deals with the /t-term of finite volume form factors and section 5 is devoted to the 
conclusions. 



2 One-particle states 

2.1 Bound-states in finite volume 

Let us consider an integrable QFT with diagonal scattering. The spectrum consists of particles 
Ai, i — 1,,..,N, with masses m, which are assumed to be strictly non-degenerate. Asymptotic 
states are denoted by 

\6 1 ,9 2 , ■ ■ • ,Qn) ili2 ... in , 

where the indices i\ . . . i n denote the particle species. Multi-particle scattering processes are de- 
scribed by the products of two-particle phase shifts Sij (#y ) where 6ij is the relative rapidity of the 
incoming particles Ai and Aj. 

The explicit formulas for the leading finite size mass corrections in a periodic box of volume L 
read [6| 

Am^) = -Y^6{ml-\ml-ml\)^ h {T c ah fe-^ L (2.1) 

AmP = -VW —e- m » Lcosh ^m bC osh(0)(S ab (e + iiT/2)-l) (2.2) 
b J ^ 2n 

\x° ah is the altitude of the mass triangle with base m c (see figure l2~!Tj) and (r^ h ) 2 is the residue of 
S a b{0) corresponding to the formation of the bound state. 
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Here we determine the leading ^-term associated to a moving one-particle state A c . Based 
on the description of mass corrections it is expected that this contribution is associated to the 
fusion A a At — > A c with the smallest \i c ah . We assume that a — b, ie. the fusion in question is a 
symmetric one. This happens to be true for the lightest particle in models with the "<I> 3 -property" 
and for other low lying states in most known models. At the end of this section we comment on 
the possible extension to nonsymmetric fusions. 

The bootstrap principle for a symmetric fusion consists of the identification 

\e) c ~\e + iui c ,e-iui c ) aa (2.3) 

resulting in 

m c = 2m a cos« c ) {fx c aa } 2 = m 2 a - -± 

Smallness of [i c aa means that m c is close to 2m a , in other words the binding energy is small. 

For a moment let us lay aside the framework of QFT and consider quantum mechanics with an 
attractive potential. Bound states are described by wave functions 

■*{x 1 ,x 2 ) = e rP ^ +x ^il>(x 1 ~x 2 ) 

where P is the total momentum and ip{x) is the appropriate solution of the Schrodinger equation 
in the relative coordinate. It is localized around x = and shows exponential decay at infinity. 
Except for the region X\ « x%, the wave function can be approximated with a product of plane 
waves with imaginary momenta p\. 2 = P ± ik. The interaction results in the quantization of the 
allowed values of k. 

The theory in finite volume is described along the same lines. There are however two differences: 

• The total momentum gets quantized. 

• ip( x ) (and therefore k) obtains finite volume corrections. 

This picture also applies to relativistic integrable theories. We consider A c as a simple quantum 
mechanical bound state of two elementary particles and use the infinite volume scattering data 
to describe the interaction between the constituents. To develop these ideas, let us consider the 
spectrum of the theory defined on a circle with circumference L. We state the identification 

\A C (9)) L ~ \Mh)Mh)) L (2-4) 

where the 0i,2 are complex to describe a bound-state; this idea also appeared in |14[ 115). Relation 
(|2.4p can be regarded as the finite volume realization of (|2.3p . The total energy and momentum of 
the bound state have to be purely real, constraining the rapidities to take the form 

01 = 6 + iu, 2 = -iu (2.5) 

where the dependence on L is suppressed. Energy and momentum are calculated as 

E = 2m a cos(m) cosh(0) p — 2m a cos(m) sinh(0) (2-6) 

The momenta of two-particle states in finite volume are quantized by a relation involving the 
scattering phase shift [TJ. This procedure can be extended in 1 + 1 dimensional integrable models 
to arbitrary multi-particle scattering states. The quantization condition for an n-particle state is 
given by the Bethe-Yang equations 

a 

e ip^Y[S ijik (e j -9 k ) = l, j = l...n 

k=l 

To quantize the bound state in finite volume, an appropriate analytic continuation of the above 
equations with n = 2 can be applied. This procedure is justified by the same reasoning that leads 
to original Bethe-Yang equations: one assumes plane waves (with imaginary momenta) except for 
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the localized interaction, which is described by the S-matrix of the infinite volume theory. Inserting 
l|2.5p and separating the real and imaginary parts 

girn a cos(u) sinh{9)L ^— m a sin(u) cosh(#)L g (2zu) 1 (2 7) 

^im a cos(u) sinh(9)L^m a sin(u) cosh(#)Lg ^ 2ili) 1 (2 8) 

Multiplying the two equations and making use of S(2iu) — S(—2iu)~ 1 one arrives at 

e 2im co S («)si n h(e)i = 1 or 2 m a cos(u) sinh(tf) = ^ (2.9) 

Li 

which is the quantization condition for the total momentum. / is to be identified with the momentum 
quantum number of A c . The quantization condition for u is found by eliminating 6 from (|2.7p : 



-m a L sin 



(«Vi+U/c'o. ( „,7 S aa (2iu) = (-1) 7 (2.10) 



The exponential factor forces u to be close to the pole of the S-matrix associated to the formation 
of the bound-state. For the case at hand it reads 

with ujj a = 2zt° c . Note the appearance of (— l) 7 on the rhs. of l|2.10p . which is a natural consequence 
of the quantization of the total momentum. This sign determines the direction from which the pole 
is approached. 

The exact solution of (|2. 10|) can be developed into a power series in e~^°- aL , where the first term 
is found by replacing u with w° c in the exponent: 

u - < c = (-l)'i (Kaf e-^vMW" + 0(e" 2 ^^) (2.12) 
First order corrections to the energy are readily evaluated to give 

E = E Q -{-lY{Yl a f^e-^ L + 0{e-^ L ) (2.13) 
where i?o is the ordinary one-particle energy 

= \ ml 



¥) 2 

In the case of zero momentum the former result simplifies to the leading term in l|2.ip . For large 
volumes we recover 

2-kI 

u — > u a nr 9 — ► arsh - 

Having established the quantization procedure we now turn to the question of momentum 
quantum numbers inside the bound state. For the phase shift let us adopt the convention introduced 
in [12] 

S ab (6) = 5„6(0)e w -W 

where S ab (6) is defined to be continuous on the real line and is antisymmetric by unitarity and real 
analyticity. 5 ab (0) can be extended unambiguously to the imaginary axis by analytic continuation 
apart from the choice of the logarithmic branch. For a generic rapidity one has 

5ab(0*) = Sab{e)* 6 ab (-6) = -S ab (0) 

Note that 5 ab (iu) is purely imaginary. 

With this choice of the phase shift the Bethe-Yang equations in their logarithmic form 

I sinh(6» + iu) + S n (2iu) = 2irl x 
I sinh(0 — iu) + Sn (— 2iu) = 2ttI 2 
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imply Ii = I 2 . Quantization of the total momentum on the other hand requires I\ = h =1/2. Note 
that different conventions for d a b would result in a less transparent rule for dividing I among the two 
constituents. The only disadvantage of our choice is the appearance of the unphysical half-integer 
quantum numbers. 

Let us denote a multi-particle state in finite volume as 

\{Il, . ■■In}) ili2 ... n>L 

where the quantum numbers Ii serve as an input to the Bethe-Yang equations. The bound-state 
quantization can be written in short-hand notation as 

\{I}) cL ^\{I/2,I/2}) aaL 

The example of mass corrections suggests that in order to obtain the total //-term it is necessary 
to include a sum in 1|2.13|1 over the different fusions leading to A c . However, the case of nonsymmetric 
fusions requires special care. Here we outline the difficulties of the nonsymmetric bound state 
quantization. 

A straightforward application of the Bethe-Yang equations to an \A a (8 + iu b lc )Ai,(9 — iu^ c )) L 

bound state yields a nontrivial phase factor e 2l7r m °-+ m b instead of (— I) 1 in <|2. 1Q[) . This in turn 
implies that the rapidity difference of the constituents can not be purely imaginary. The real parts 
of the rapidities thus get different finite size corrections and the total energy of the bound state 
becomes complex. A possible solution would be to take twice the real part of the energy correction, 
corresponding to the sum of the contributions coming from the A a A b — ► A c and A b A a — > A c fusions. 
This is however only a guess and a more systematic treatment is needed. In fact, the infinite volume 
bootstrap principle suggests that the states 

\A a (9-iu b ac )A b (e + iuZ c )) L and \A a (9 + iu b ac )A b (9 - m£ c )) £ 

should be handled on an equal footing. A possible way to accomplish this would be to develop a 
multi-channel Bethe-Yang quantization scheme. However, this is beyond the scope of the present 
work. 

As a conclusion of this section (|2- 12|) is compared to the lowest order results of the TBA 
approach. The general discussion of excited states TBA equations in diagonal scattering theories 
is not available. For simplicity we restrict ourselves to the Lee- Yang model which was considered 
in the original paper [16]. In this nonunitary model there is only one particle and the scattering is 
described by 

= sinh(fl) +zsin(7r/3) 
1 ' sinh(0) -isin(7r/3) 
The exact TBA equations for moving one-particle states read 

/OO 
— mcosh(0)Z(0) (2.14) 
-oo 27r 

e(9) = mflcoshfl + log S( f ~ 9o \ -(if* L){6) (2.15) 

/OO 1/1/ 
— f(9-9')g(8') 



where 



Here the volume is denoted by R to avoid confusion with L(9). The complex rapidity 9q satisfies 
the consistency equation 

e(9 ) = mR cosh 9 + in - \og(S(2Hm9 )) - {<p-kL)(9) = i{2n + 1)tt (2.16) 

The convolution term in (|2 . 16[) can be neglected and one obtains Im#o = tt/6 + 5 where S is 
exponentially small. To zeroth order one also has 

mR cosh (Re6> ) = (in + (1 — sign<5))7r 
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which is the ordinary one-particle quantization condition with I = 2?i + i(l — sign<5). Neglecting the 
contribution of the integral in <|2. 14|) and substituting 6 = 9 + iu one has E — 2m sin(w) cosh(#) . 
This is exactly the energy of an AA bound state with the imaginary rapidities 8±iu. Separating the 
real and imaginary parts of (|2.16j) and still neglecting the convolution term (which is responsible 
for the F-term) one obtains equations (|2 . 10[) and (|2,9j) . thus proving the consistency of the two 
approaches. 

It seems plausible that by exactly solving the bound state quantization condition (|2. 10|) one 
obtains all higher order corrections that go as e~ nf * aaL with n e N. In the next subsection we 
present numerical evidence to support this claim. 

2.2 Numerical analysis 

We investigate the famous E$ scattering theory [17], which is the relativistic integrable field 
theory associated to the scaling limit of the Ising model in the presence of a magnetic field. The 
infinite volume spectrum of the model consists of 8 particles. The first three particles lie below the 
two-particle threshold and they all show up as AiAi bound states. These fusions are responsible 
for the leading //-term. The corresponding parameters (in units of mi) are listed in the table below. 
The exponent of the next-to-leading correction (the error exponent) is denoted by //. 



a 


m a 


a 
Mil 






1 


1 


0.86603 


205.14 


1 


(mi) 


2 


1.6180 


0.58779 


120.80 


0.95106 


(M?a) 


3 


1.9890 


0.10453 


1.0819 


0.20906 


(2Mii) 



In |6J Klassen and Melzer performed the numerical analysis of mass corrections. The analytic 
predictions were compared to TCSA data and to transfer matrix results. They observed the expected 
behaviour of mass corrections of Ai and A% in the former case they were also able to verify the 
F-term. On the other hand, the precision of their TCSA data was not sufficient to reach volumes 
where the /j,-term for A3 could have been tested. This limitation is a natural consequence of the 
unusually small exponent /Ltf j : the next-to-leading contribution is of order e~ 2 ^ llL , still very slowly 
decaying. 

Here we employ the TFCSA (Truncated Fermionic Conformal Space Approach, see [18]) routines 
that were successfully used in [9J [12] . Calculations are performed for J = 0, 1, 2, 3 at different values 
of the volume. One-particle states of A±, A2 and A3 are easily identified: they are the lowest 
lying levels in the spectrum, except for 1 = where the lowest state is the vacuum. We use the 
dimensionless quantities I = m\L and e = Ejm\. 

The results are extrapolated from e cut = 20.. 30 to e cu t = 00 using the procedure developed 
in [19]. Our experience shows that this extrapolation technique reduces the numerical errors by 
an order of magnitude. However, our attempts to develop an adequate method to estimate these 
errors have failed. Several examples were encountered where the actual numerical deviation was 
either underestimated or overestimated, no matter which estimate was used. We therefore resign 
from quantitatively monitoring the TCSA errors and constrain ourselves to a range of the volume 
parameter where it is safe to neglect truncation effects. 

2.2.1 A 3 

We begin our analysis with the most interesting case of A3 . At each value of I and / the following 
procedure is performed. 

• The energy correction is calculated according to (|2.13fl 

• The quantization condition (|2 . 10[) is solved for u and the energy correction is calculated by 

Ae = 2 cosh((9) cos(it) — eo 

where is determined by the total momentum quantization (|2.9p and e is the ordinary 
one-particle energy. 

• The exact correction is calculated numerically by Ae = e TCSA — e . 
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Figure 2.2: Finite size corrections to A3 one-particle levels (in sectors 
I = ... 3) as a function of the volume. The TCSA data are plotted against 
theoretical predictions of the single /i-term associated to the A\A± — > A3 
fusion (solid curve) and the exact solution of the bound-state quantization 
(dotted curve). 



The choice for the range of the volumes is limited in two ways. On one hand, I has to be 
sufficiently large in order to reduce the contribution of the F terms and other higher order finite 
size corrections. On the other hand, numerical errors grow with the volume and eventually become 
comparable with the finite size corrections, resulting in an upper bound on I. The window I — 30.. 40 
is suitable for our purposes. 

The results are shown in figure [221 It is clear that the /x-term yields the correct prediction in 
the L — > 00 limit. However, higher order terms cause a significant deviation for I < 40, which is in 
turn accurately desribed by the bound state prediction. The sign of the correction depends on the 
parity of I as predicted by (|2 . 13[) . 

Based on the success of this first numerical test we also explored the region I < 30. Inspecting the 
behaviour of u as a function of I reveals an interesting phenomenon. It is obvious from (|2. 12|) that 
u(l) is monotonously increasing if I is odd, with the infinite volume limit fixed to u a ac . However, the 
complex conjugate pair 6*1,2 approaches the real axis as I is decreased and they collide at a critical 
volume I — l c . For I < l c they separate again but stay on the real line, providing a unique solution 
with two distinct purely real rapidities. The same behaviour was also observed in [14} 115). 

The interpretation of this phenomenon is evident: if the volume is comparable to the character- 
istic size of the bound-state, there is enough energy in the system for the constituents to become 
unbound. Therefore the A3 one-particle level becomes an A\A\ scattering state for I < l c . We call 
this phenomenon the "dissociation of the bound state". The same result was obtained also in the 
boundary sine-Gordon model by a semiclassical analysis |20j . 

The value of l c can be found by exploiting the fact that the Jacobian of the Bethe-Yang equations 
(viewed as a mapping from (61,62) to (h^h)) vanishes at the critical point. A straightforward 
calculation yields 

i c = VWO) 2 -/ 2 ^ 2 

where fu(6) = 8^(9). The numerical values for the case at hand are 

l c = 27.887 (I = 1) and l c = 26.434 (/ = 3) 
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Figure 2.3: A 3 one-particle levels in sectors / = ... 3 as a function of the 
volume. Dots represent TCSA data, while the lines show the corresponding 
prediction of the AiAi bound state quantization. In sectors I — 1 and 
1 = 3 the bound state dissociates at l c and for I < l c a conventional A\Ai 
scattering state replaces A3 in the spectrum. The values of l c are shown 
by the two arrows. 

We are now in the position to complete the numerical analysis. The Bethe-equations are solved 
at each value of I, providing two distinct real rapidities for / < l c (with / being odd), and a complex 
conjugate pair otherwise. The energy is calculated in either case as 

e = cosh(#i) + cosh(^2) 

which is compared to TCSA data. The results are exhibited in figure 1231 

The agreement for the upper two curves (1=1 and / = 3) is not as suprising as it may seem 
because what one sees here are conventional A1A1 scattering states. The Bethe-Yang equation 
determining their energy is exact up to 0(e~^ L ) where fi' = ls the smallest exponent that 
occurs in the sequence of finite volume corrections of A\. On the other hand, the energy levels are 
analytic functions of L, which leads to the conclusion that the prediction of (|2. 10|) is correct up to 
even for L > L c . Comparing the numerical values one finds > 8/1^. We conclude 
that the bound state picture indeed accounts for finite volume corrections up to the first few orders 
in e _M <"' i (the first 8 orders in the case at hand). 

2.2.2 A x and A 2 

Particles Ai and Ai also appear as A\A\ bound states. However, there is no point in applying 
the complete bound state quantization to them, because the error terms dominate over the higher 
order contributions from <|2. 1Q[) : the exponents of the subleading finite size corrections mi and 
H12 are smaller than and 2/i.f x . Nevertheless, the leading /i-term can be verified by choosing 
suitable windows in I. 

In figures HOI and [A.2I logf| Ae|) is plotted against the prediction of 1|2.13|1 for I = 6. .16 and 
I = 6. .22. (the sign of Ae was found to be in accordance with (|2- 13[) for both A\ and A2) 

In the case of A\ perfect agreement is observed for I = 10. .18 in the sectors 1 = and 1=1. 
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For 1 = 2 and 1 = 3 the energy corrections become too small and therefore inaccessible to TCSA 
(note that the prediction for I = 3 is of order 10~ 6 ). 

In the case of A 2 precise agreement is found for I = 14. .22 in all four sectors. 

2.2.3 A 5 

Here we present an interesting calculation that determines the leading mass corrections of 
A$. The standard formulas are inapplicable in this case, because lies above the two-particle 
threshold. However, it is instructive to consider the composition of A§ under the bootstrap principle 
and to evaluate the /x-term prediction. 

There are two relevant fusions 

A 1 A 3 , -> A 5 with i4 3 = 0.2079 
A 2 A 2 -> A 5 with [i\ 2 = 0.6581 

Numerical evaluation of (|2.ip shows that the contribution of the second fusion is negligible for 
/ > 30. The first fusion on the other hand yields a significant discrepancy when compared to TCSA 
data. This failure is connected to the two-particle threshold and it can be explained in terms of the 
bound state quantization. Experience with A3 suggests that one should first take into account the 
energy corrections of ^3 and consider the A1A3 — > A5 fusion afterwards. ^3 can be split into A\A\ 
leading to the "triple bound state" A\A\A\ — > A5. In infinite volume one has (see also fig. 12.11 c.^l 

|0) B ~ |fl-2»u? 1 ,fl,fl + 2*fi? 1 > 111 

The finite volume realization of this identification is most easily carried out in the 1 = sector 
with 

|{0}> 6l z~K0,0,0}> llliZ 
Setting up the three-particle Bethe-Yang equations with rapidites (iu,0, —iu): 

e- miSin ^ L S n (iu)S n (2iu) = 1 
Sii(iu)Sn(—iu) = 1 
e misin(u)i Sii(-m)Sii(-2m) = 1 

The second equation is automatically satisfied due to unitarity and real analyticity, whereas the 
first and the third are equivalent and they serve as a quantization condition for u. The finite volume 
mass of A§ is given in terms of the solution by 

m 5 (l) = 2cos(u) + 1 (2.17) 

In the large L limit the infinite volume mass is reproduced by u — > 2ttf 1 . Figure l2~4l demonstrates 
the agreement between TCSA and the prediction of H2.17J1 . 

The possibility of solving the quantization of the triple bound state in a moving frame looks 
very appealing. In the general case the rapidities are expected to take the form (9\ +iu, 62, 61 — iu) 
where 9\ and 9 2 do not necessarily coincide. However, the numerical precision of our TCSA data 
was not sufficient to check our predictions. 

3 Multi-particle states 

3.1 Bethe-Yang quantization in the bound state picture 

Multi-particle states in finite volume are denoted by 

where the quantum numbers Ij serve as an input to the Bethe-Yang equations 

...,9 n ) = m l3 sinh(0j)£ + ^ S ijih (9j - 9 k ) = 2-nIj j = l...n (3.1) 
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Figure 2.4: Finite size mass of A§ as a function of the volume. The 
squares represent the TCSA data which is compared to the leading \i- 
term (solid curve) and the solution of the quantization condition for the 
"triple bound state" ^4i^4iAi (dotted curve). The straight line shows the 
infinite volume mass. 



The energy is calculated as 

n 

E = ^2m ij coaix(dj) + ... (3.2) 
3=1 

where (S, Si,... ,S n ) is the solution of l|3.ip . The dots indicate exponentially decaying finite size 
corrections. 

Let us consider a scattering state \{I, h,. ■■ , In}) cbl & l composed of n+ 1 particles, the first 
one being A c . We determine the leading part of the [i-term by considering A c as an A a A a bound 
state inside the multi-particle state. Therefore we write 

\{I, In}) cblM ~ |{//2, 1/2, h,..., I n }) aablM (3.3) 

The energy is then determined by analytic continuation of the n + 2 particle Bethe-Yang equations. 
They read 

n 

e -^aSU,(u)cosh(9)L etma cos(u)sinh(9)L Sa ^ 2lu ^S ab] (e + lU-6 :j ) = 1 (3.4) 

3=1 
n 

e ™*sMu)co S h(9)L e t ma cos(u) S inh(e)L Saa (_ 2iu jY[S ab] (8~ iu - 6 j) = 1 (3.5) 



3 = 1 



e im bj sinh(0j)L 



Sab^ej-e-i^Sab^ej-e + i^USb^iej-ek) = 1 (3.6) 



k=l 



The ordinary n + 1 particle Bethe-equations are reproduced in the L — > oo limit by multiplying 
l|3.4p and l|3.5p and making use of the bootstrap equation 

S c bM=S abj (9 + iul c )S abj {e-iu a ac ) 

We now proceed similar to the previous section and derive a formula for the leading correction. 
The shift in the imaginary part of the rapidity can be calculated by making use of l|3.4p and (|2,lip 

as 

(yc \ 2 _ n 

All = U-U a ac = L^L e -Mcosh(e)L elmc si„h(0)L/2 "Q (g + .-^ _ ^ ^ 

3=1 
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Multiplying $3A§ and flU} 

n 

e l2m " cos(u) sinh(e)L J] 5 afcj (6-iu- 9 )S abj {e + iu-Oj) = 1 



^im^. sinh(#j)L 



S«6,(fli-e-iu)S , i i (fl i -fl + iu)ij56 i 6 lk (ei-dfc) = 1 



fe=l 



Let us define 
where 



S abj (6-iu- 9j)S abj {9 + iu- dj) » 5 chj (0 - 9j)e iAu ^ b i 



^ cfcj (#) = i<p cbj [0 + iu a ac ) - itp cbj (9 - iu a ac ) with <p ab {9) = 6' ab {9) 
Using 2m a cos(w) w m c — 2^i c aa Au the logarythm of (|3.8p and (|3.9p can be written as 



(3.8) 
(3.9) 



Q j (9,9 1 ,...,9 n ) = $ cbj {9-9 )Au 
The lhs. can be expanded around the n + 1 particle solution (9,9 n , . . . , 9 n ) to arrive at 



/ - \ 

0i -Si 

\#n — 9 n / 



J 



(n+1) 



/2 M S„ sinh(0)£ - £?=i £ c6 . (0 - 0,)\ 



\ 



<p cbn (9 - 9 n ) 



Ai 



J 



where 



The final result for the energy correction reads 



AE = -2/4, cosh(0)Aw + 



/ m c sinh(0) \ 
m bl sinh(0i) 



\m bn sinh(0 n )/ 



(3.10) 



(2^ aa sinh(0)£ - E"=i <Po bj (9 - 9 3 )\ 
£ cbl (9 -h) 



V 



fcb n {9 ~ 9 n ) 



A; 



(3.11) 

with Au given by (|3~7j) . 

Based on the previous section it is expected that there is a similar contribution for every fusion 
leading to each one of the constituents of the multi-particle state. 



3.2 Multi-particle states - Numerical analysis 

We first consider finite size corrections to A\ A3 states. They are not the lowest lying two-particle 
states in the spectrum, but they possess the largest /i-term which is connected to the A1A1 — > A3 
fusion. Given a particular state \{Ii,h})i3 L the following procedure is performed at each value of 
the volume: 

• The two-particle Bethe-Yang equation for \{Ii, is}) 13 L is solved and the /i-term is calculated 
according to (|3.1ip . 

• The exact three-particle Bethe-Yang equation is solved for |{/i, I3/2, /3/2}) m L 

The results for different Ax A3 levels are shown in figure HOI The situation is similar to the case 
of the A3 one-particle levels: the bound state quantization yields a remarkably accurate prediction, 
whereas the single /i-term prediction only becomes correct in the L — * 00 limit. 
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In table IA.ll we present a numerical example for the dissociation of the bound state inside the 
two-particle state. In this case an AiA 3 state turns into a conventional A1A1A1 three-particle state 
at l c w 30. 

Finite size corrections to A\A± and A1A2 states are also investigated, the leading /i-term given 
by the fusions A\A\ — > A± and A\A\ — > A2, respectively. In the former case we calculate separately 
the contribution associated to both A\ particles and add them to get the total correction. Results 
are exhibited in figures lA~4l and lA~5l and formula (|3.1ip is verified in both cases. 

4 Finite volume form factors 

The connection between finite volume and infinite volume form factors was derived in [I2J as 

h...j m ,L{{l'v • • • , OP(0, 0)|{/l, . . . , M> n ..,„,L = 

+ O(e -/1 ' i (4.1) 



F°(9' m + iTr,...,6' 1 +in,9 1 ,..., 9 n ) ]m ... hll ... ln n(p -^L\ 



Pn...i n (dl, ■ ■ ■ ,^n)Pj 1 ...j m (d'i, ■ ■ ■ ,8'rn) 

where the rapidities 9 are solutions of the corresponding Bethe-Yang equations and it is supposed 
that 9j ^ 9' k whenever ij = ik- The extension of l|4.ip to include disconnected terms can be found 
in [9] . The proportionality factor in (|4,ip is given by the density of states 

pM in (§ 1 ,...,6 n ) = detjW , J^ = ^ , k,l = l...n 

which can be used to identify formally the finite volume and infinite volume states as 

\{h, ■ ■ ■ In}) ili2 ... n x ===^===\§ 1 , dn) n _ ln 

Based on general arguments it was shown in [12] that [i' > pi where /i is determined by the pole of 
the S-matrix closest to the physical line. A systematic finite volume perturbation theory (Liischer's 
method applied to form factors) is not available. However, it is expected that the actual value of 
fi' depends on what diagrams contribute to the form factor in question. Apart from the insertion 
of the local operator they coincide with the diagrams determining the finite size corrections of the 
multi-particle state. Therefore p! is associated to the bound state structure of the constituents of 
the multi-particle state. In this section it is shown that the leading correction term can be obtained 
by the bound state quantization. 

4.1 Elementary one-particle form factors 

(|4.ip yields a simple prediction for the elementary one-particle form factor: 

Ff (I, L) = (0\O(0, 0)\{I}) c r =^L + 0(e-»' L ) (4.2) 



where E is the one-particle energy, and F® = F®{9) is the infinite volume one-particle form factor, 
which is constant by Lorentz symmetry. 

The p-term associated to (|4.2p is derived by employing the bound state quantization. We gain 
some intuition from the previous sections where it was found that the bound state A a A a may 
dissociate at a critical volume L c . For L < L c there is no one-particle level of type A c in the given 
sector of the spectrum, however an A a A a scattering state appears instead. Finite volume form 
factors of this state are calculated using (|4.ip as 

Ff(/,L) = ^%% forL<L c (4.3) 

The generalization to L > L c seems to be straightforward: one has to continue analytically (|4.3p 
to the solutions of the Bethe-Yang equation with imaginary rapidities #1,2 — 9±iu. However, note 
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that equations l|4.2p and l|4.3p are valid up to a phase factor. In order to continue analytically to 
imaginary rapidities we also need to fix this phas^]. 
The two-particle form factor satisfies 

F°(9i, 9 2 )aa = S aa (9i — 92)F°(02, 0\) aa 

The simplest choice for the phase is therefore 



F°(8 1 ,8 2 ) aa = ^S aa {9 l -e 2 ) \F°(d 1 ,e 2 ) aa \ (4.4) 



This choice is dictated by CPT symmetry [21j [T0|. fTT] , and it is respected by all known solutions of 
the form factor bootstrap axioms. There is a sign ambiguity caused by the square root, but it can 
be fixed by demanding (S' aa (0)) 1 / 2 = i and continuity. Using l|4.4p 



-<0(T T\ — V $00(62 — 6\)F (61,62)00 



F?(I,L) = 



and upon analytic continuation 



^ r ^_ ^S aa (-2iu)F°(9 + iu,6- iu) aa 



F?(I,L) = v "" v ' \ ' '— for L>L C (4.5) 

y/p aa (0 + iu,9 — iu) 

It is easy to see that the result l|4.2p is reproduced in the L — > oo limit. First observe that 

F°(6 + iU,6- iu) aa ~ Vla _ r r Jf(g) 
2(U ~ Ka) 

The residue of p aa is determined by ip aa (2iu) and it reads 

S' (2m) 1 
p{6 + iu,6 - iu) aa ~ 2m a L cos(u) cosh(#) (-£) aa ) = m c L cosh(#) ; 



S aa (2iu) ' ' 2(« — SSJ 

The singularities in the numerator and denominator of l|4.5p cancel and indeed 

F?{I,L) 



\Jm c L cosh(#) 



We emphasize that it is crucial to include the extra normalisation factor \J S aa (— 2iu) to obtain a 
meaningful result. 

Expression l|4.5p can be developed into a Taylor-series in u — u^ c . The first order correction is 
evaluated in Appendix A and it reads 

If (I, L) = -?L= - 2iF "°(g) (u _ 6 » ) + 



F? 

+- 



2££° , ™cM f m l E _ | L 



(rSa ) 2 (E°r \m* c E° 



(u - u a ac ) + 0(e-^ L ) (4.6) 



where 



S C J = lim (S qq (2j U ) "" J 



2(« - t*Sc) / 

is the ordinary one-particle energy and the rapidity shift u — u a ac is given by l|2.12p . 



x The phase of a (nondiagonal) infinite volume form factor is unphysical in the sense that it may be redefined by a com- 
plex rotation of the state vectors and physical quantities, e.g. correlation functions, do not depend on such redefinitions. 
However, the bootstrap program uniquely assigns a phase to each form factor. 
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4.2 One-particle form factors - Numerical analysis 

Let us introduce the dimensionless form factors as 

")\e(0,0)\{I}). L 



Mi, i) 



mi 



The machinery of [12] is used to determine /*(/, I) for i = 1,2,3 and I = 0, 1, 2, 3. The numerical 
results are compared to the exact infinite volume form factors [22j [23] . 

We start our investigation with fs(I,l), for which relatively large exponential corrections were 
already reported in [12]. It is convenient to consider 



h{I,l) = {e l) 1/2 h{I,L) with lim f 3 (I,l)=F 3 (4.7) 

l — KX) 

The numerical results are demonstrated in fig. 14.11 Note, that this is exactly the same figure as 
4.6. (c) in [12J, but this time the interpretation of the huge deviations from F 3 is also provided. 

We also tried to verify the predictions for /1 and fa . In the latter case reasonably good agreement 
was found with TCSA, the results are demonstrated in fig. IA.6I In the case of /1 we encountered 
the unpleasant situation that the F-term decays slower than the TCSA errors grow, thus making 
the observation of the /z-term impossible. 

It is straightforward to generalize l|4.5p to matrix elements between two different one-particle 
states. For b =^ c one has for example 

/rj -,, , rj - lx F e (9 b + iir,9 + iu,9-iu) baa 



pb{0b)Paa(0 + iu,6— iu) 

Numerical examples are presented in figures [4721 (a)-(c) for c = 3 and b = 1,2. 

The most interesting case is the one shown in fig. 14.21 (d) where the matrix element between 
two different A3 one-particle states are investigated. This can be done by considering both A3 
particles as the appropriate A1A1 bound states and then calculating the finite volume form factor 
3<{J}|e|{r}> s , £ as 

n ({J/2, 1/2}\e\{l'/2, l'/2}) n , = ^ + ^ + ^,9 -iu + i^ff +iu',9' -iu'^u 

u 1 "".J y/ Pll {e> + iu>, 9' - iu') Pll {o + iu,e- iu) 

Once again we find complete agreement with the TCSA data. 

4.3 Elementary multi-particle form factors 

The generalization of l|4.5p to multi-particle states is straightforward, the only task is to find 
the appropriate phase factor. Similar to the one-particle case one has 



F°(e 1 , . . . , e m ) h ... bm = /JJ S bibj (9, - 9j) \F°(9 X , . 



?m )bi ...6„ 



A general n particle finite volume form factor with real rapidities can thus be written as 



/ IL <3 S btb] (9 3 - 0j) 
y p n (9i, . . ■,9 n )b 1 ...b„ 

Substituting the solution of the Bethe-equation for the state \{I/2, 1/2, h, . . . , In}) aabl b L and 
making use of the real analycity condition 

\S abj (0j - 9 - iu)S ab] (9j -6 + iu)\ = l 

one gets 

{0\U\{l/2,l/2,h,...,l n }) aab b L / mmfXT ■ R ■ a a \ { ' 

V 'p( n + 2 ) (9 + IU,9~ IU, 61,..., 9 n )aab 1 ...b n 
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/ 
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35 
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Figure 4.1: Elementary finite volume form factors of A3 one-particle lev- 
els. Here the normalization (|4,7p is applied to obtain a finite I — » 00 limit, 
which is given by the infinite volume form factor F3 = (0|e|^43(6')). The 
TCSA data are plotted against the bound state prediction. The ordinary 
evaluation of fa is simply the constant F 3 . 




Figure 4.2: One-particle-one-particle form factors, dots correspond to 
TCSA data. The solid lines represent the ordinary evaluation of the fi- 
nite volume form factors, while the dotted lines show the bound state 
prediction. 
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0.24 





20 22 24 26 28 30 32 

(b) <0|e|{2,0}> 13!i 



Figure 4.3: Elementary form factors of AiA 3 scattering states, dots 
correspond to TCSA data. The solid lines are obtained by a "naive" eval- 
uation of the finite volume form factors, while the dotted line represents 
the bound state prediction, (in this case A1A1A1 form factors at the ap- 
propriate rapidities) 



up to a physically irrelevant phase. 

It is easy to show once again that the "naive" result is reproduced in the L — > oo limit. To do 
so, we first quote the dynamical pole equation of the infinite volume form factor: 

F°(6 + m, - iu, ft, ... , 0„) oabl ... 6 „ = F ° Q _ a , F°(9, 0i, ... , 6 n ) cbl ... bn + 0(1) 

2(u - Kc) 

The singularity of p( n+2 ) is given by 

Res u ^ <c(O (n+2) (0 + iu, e - iu, 0i, ... , e n ) aabl ... bn = ^ (n+1) (0, 0i, • • ■ , e n ) cbl ... bn 

The "naive" formula is now recovered by inserting the last two equations into l|4.8p . 

The leading exponential corrections can be obtained by plugging (|3.7p and p.lOp into (|4.8[) and 
expanding to first order in u — u a ac . This procedure is straightforward but quite lengthy, therefore 
we refrain from giving the details of the calculations. 

In fig. 14.31 two examples are presented for the evaluation of 1|4.8|1 applied to A\A 3 two-particle 
states. 



5 Conclusions 

In this work we determined the generalization of Liischer's /i-term associated to moving one- 
particle states, arbitrary scattering states and finite volume form factors. Our method is based on 
the bootstrap principle of the infinite volume theory which states that a particle which is a bound 
state of two others is indeed indistinguishable from the two-particle state with the appropriate 
(imaginary) momenta. An analytic continuation of the Bethe-Yang equations was used to quantize 
the bound states in finite volume. 

The analytic results were tested by comparing the predictions to high-precision TCSA data of 
the Ising model; a satisfactory agreement was observed in each case. We also demonstrated that the 
bound state quantization goes beyond the leading mu-term corresponding to the fusion and gives 
a resummation of all powers of e-i*™ L . This is a substantial improvement over the inclusion of the 
leading Liischer term when other sources of finite volume corrections (other fusions, F-terms) can 
be neglected. Our main example was provided by particle ^3 from the Es scattering theory. 

As an analytic check the calculations were compared to the leading order results of the exact 
TBA equations. The agreement between the two approaches was explicitly demonstrated in the 
case of the one-particle levels of the Lee- Yang model. 
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Our results can also be applied in nonintegrable models for states below the first inelastic 
threshold. However, the calculations only apply to symmetric fusions of the type A a A a —> A c . 
The nonsymmetric case requires an extension of the Bethe-Yang equations to incorporate multi- 
channel scattering. Such a quantization scheme could also be used to describe finite volume states 
in nonintegrable theories above the two-particle threshold. 
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A ^-term for the one-particle form factor 

Here we develop the first order correction to (|4,5|1 . Using the exchange axiom 

F°{B-iu,6 + iu) aa 



F C U (I,L) 

The form factor axioms imply that 

lim F°( 



y/ S aa {-2iu)p aa (B + iu,9~ iu) 



(A.1) 



-F. 



The simple pole of <p aa (2iu) in p aa is cancelled by S aa (—2iu), therefore both the numerator and 
the denominator of (|A,1|) have continuous limits asw-> SJ C . 

The form factor F°{9 — iu,9 + iu) aa only depends on u by Lorentz-symmetry. Therefore 

F°{9 -iu,6 + iu) aa = ±-F? - 2i (Fg)' (u - u a ac ) + ... 



where 



0-e r =-2iu% a 

Expanding the S-matrix element into a Laurent-series in the vicinity of the pole 



S aa (2iu) = 



(ru 2 
2(« - Ko) 


, cc,0 


2(« - Ko) 







S aa (-2iu) 

Expanding the denominator: 

S(—2iu)p aa (9 + iu,9 — iu) 



OC,0 , 
"».n I 



S{-2iu)E l E 2 L' 1 - i{E x + E 2 )LS'(2iu)[S(-2iu) 



E C L 2 El E 2 L* 4E C L 



where 



Putting all this together 



E c = Ei + E 2 = 2m a cos(w) cosh(#) 



F° 



Ve~l 



-2^ (FZ) 



E 



x/ej: 3 



-EiE 2 L 2 + — ° s2 Saa 



(rs a ) 
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Note that in the preceding formulas E c does include the leading order correction to the usual 
one-particle energy E° = v / m 2 + (2nI) 2 /L 2 . Using 

E c = E° c - 2^(« - Kc) + 0{e-^ L ) and E,E 2 = ^E 2 fj 2 



rri 



the final result is given by 
F C °(I,L) = 



y/E°L yfEfL 



(« - <c) + 



+ ■ 



H§ rO - £ ) L 



<c) + 0{e- 2 » L ) (A.2) 



with 
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(a) I = 
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(b) / = 1 



10 12 14 



(c) I = 2 



(d) 7 = 3 



Figure A.l: Finite size corrections to A\ one-particle levels in sectors 
I = ... 3, logwAe is plotted as a function of the volume. Dots represent 
TCSA data, while the lines show the ^(-term corresponding to the A\A\ — > 
fusion. 



(a) 7 = (b) 7 = 1 



(c) 7 = 2 (d) 7 = 3 

Figure A. 2: Finite size corrections to A2 one-particle levels in sectors 
I = ... 3, logwAe is plotted as a function of the volume. Dots represent 
TCSA data, while the lines show the ^-term corresponding to the A X A X — » 
A 2 fusion. 
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Figure A. 3: Finite size corrections to A1A3 scattering states as a func- 
tion of the volume. Dots represent TCSA data, the solid line shows the 
/x-term corresponding to the A\Ai — > A3 fusion. The dotted lines are ob- 
tained by the exact solution of the quantization condition for the A1A1A1 
three-particle system. 



(a) K-0.5,1.5}), 



(b) |{-0.5,1.5}) n 





(c) |{0.5,1.5}> 1 



(d) |{0.5,2.5}> n 



Figure A. 4: Finite size corrections to A 1 A 1 scattering states, log w Ae is 
plotted as a function of the volume. Dots represent TCSA data, while 
the solid line show the sum of the two /z-terms corresponding to the 
A1A1 — > Ai fusions. 
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(a) |{0,1}) 12 



GO 



Figure A. 5: Finite size corrections to A\Ai scattering states, logioAe 
is plotted as a function of the volume. Dots represent TCSA data, while 
the solid line show the /z-term corresponding to the A\A\ — > A 2 fusion. 



I 


9i 


9 2 


3 


e(l) (predicted) 


e(i) (TCSA) 


22 


0.44191 


0.68235 


-0.58742 


3.51877 


3.51900 


23 


0.42574 


0.64280 


-0.55190 


3.46201 


3.46219 


24 


0.60523 


0.41155 


-0.51892 


3.41239 


3.41255 


25 


0.56934 


0.39931 


-0.48831 


3.36890 


3.36907 


26 


0.38910 


0.53475 


-0.45989 


3.33071 


3.33089 


27 


0.50092 


0.38119 


-0.43350 


3.29709 


3.29729 


28 


0.46686 


0.37634 


-0.40899 


3.26743 


3.26767 


29 


0.42947 


0.37739 


-0.38622 


3.24122 


3.24164 


30 


0.38646 + 0.02332 i 


0.38646 - 0.02332 i 


-0.36505 


3.21801 


3.21832 


31 


0.37059 + 0.04042 i 


0.37059 - 0.04042 i 


-0.34537 


3.19741 


3.19775 


32 


0.35572 + 0.05104 i 


0.35575 - 0.05104 i 


-0.32706 


3.17908 


3.17945 


33 


0.34182 + 0.05891 i 


0.34182 - 0.05891 i 


-0.31002 


3.16275 


3.16313 


34 


0.32876 + 0.06513 i 


0.32876 - 0.06513 i 


-0.29415 


3.14816 


3.14856 


35 


0.31650 + 0.07022 i 


0.31650 - 0.07022 i 


-0.27936 


3.13511 


3.13547 


36 


0.30498 + 0.07446 i 


0.30498 - 0.07446 i 


-0.26556 


3.12341 


3.12235 



Table A.l: An example for the dissociation of the A1A1 bound state 
inside a scattering state. |{2, 0}) 31 L is identified with |{1, 1, 0}) m L and 
the corresponding Bethe-Yang equations is solved. For I < 30 there is a 
real A1A1A1 three-particle state in the spectrum, whereas at I w 30 two of 
the rapidities become complex and the two-particle state A1A3 emerges. 
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Finite size corrections to the elementary form factors of A2 . 



Dots represent TCSA data, while the lines show the ^t-term prediction 
corresponding to the A\A\ — > A2 fusion. 
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